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We investigate the collision of two oblique dark solitons in the two-dimensional supersonic non- 
linear Schrodinger flow past two impenetrable obstacles. We numerically show that this collision 
is very similar to the dark solitons collision in the one dimensional case. We observe that it is 
practically elastic and we measure the shifts of the solitons positions after their interaction. 



. 1. Two-dimensional (2D) oblique dark solitons are unstable with respect to transverse perturbations jl|-l3| and 
' therefore their interaction with each other is not of much interest from practical point of view. However, it has been 
found [3] that such solitons generated in the flow of an atomic Bose-Einstein condensate (BEG) past an obstacle 
behave as effectively stable. Such a behavior was explained in [5i] as a result of the transition from absolute instability 
J> ] of 2D solitons to their convective instability for large enough velocities of the flow in the reference frame attached 
O i: to the obstacle, so that unstable modes are convected by the flow along the solitons from the region around the 
obstacle. The condition for convective instability of these dark solitons was derived in Q. This phenomenon has a 
general nature and its nonlinear optics counterpart has been discussed in Q- Recently, experimental observations in 
Bose-Einstein condensate of exciton-polaritons have indeed demonstrated the existence of stable oblique dark solitons 
in a superfluid flow past an obstacle Q. Hence, interaction of such effectively stable oblique dark solitons becomes a 
question of considerable interest and it will be addressed in this paper. 
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^ 2. Oblique dark solitons in a superfluid are described very well |^ as stationary solutions of the defocusing nonlinear 

Schrodinger equation (NLS) 

^ ■ z^t^-iAV. + lV'PV, (1) 

>'.... 

, which is written here in standard dimensionless units and A = 9^ -I- 9^. Its transformation to a "hydrodynamic form" 
(~~«. ' by means of the substitution 

cn 
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7/;(r,t) = v/n(r,i)exp \^ J u{r',t)dr'j (2) 
yields the system 

nt + V ■ (nu) = 0, (3) 



where n is the density of the fluid and u denotes its velocity field. 
In a stationary case (n^ = 0, Ut — 0) this system takes the form 



uu, + vuy + n,+ \ - 0, 



{nu)x + {nv)y = 0, 
-0, 

(5) 
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FIG. 1: Interaction of two oblique dark solitons generated in the flow of a superfluid past two impenetrable obstacles. The flow 
is from the left to the right with M — 5. One obstacle is located at (a;io, t/io) = (0, —10) and the other one at (1:20, J/20) = (0, 10) . 
Colors indicate the magnitude of the density. 



where we have introduced the components u — (u, v) of the velocity field. It should be solved with the boundary 
conditions 

n = 1, u — M, V — at |x| — >■ 00, (6) 

which means that there is a uniform flow of a superfluid with constant velocity u — (M, 0) at infinity. Since in our 
dimensionless units the sound velocity at infinity is equal to unity, the incoming velocity coincides with the Mach 
number M . The soliton solution of this problem was found in 4] and it can be written as 

. , ^ l-M^sin^e 

n[x,yj = 1 ^ , (7j 

cosh^ [ V 1 - sin^ 9{x sin - y cos 9 + yo cos 6)] 



u{x, y) = M 



1 + sin^ 



n{x,y) 



v{x,y) 



-M sin 6* cos 6' 



n{x,y) 



(8) 



where 9 is the angle between the oblique soliton and the horizontal axis and yo is its intersection point with the 
y axis. The transformation ([2]) implies that the flow is potential (vorticity free) so that the velocity field u can be 
represented as a gradient of the phase 



(a;, y) = Mx — arctan ■ 



M sine* 



Vl — sin^ 9 tanh[-\/ 1 — fv'P sin^ 9{x sin9 — y cos 9 + yo cos 9)] 
Correspondingly, the wave function of the oblique soliton reads 

y) = I Vl - sin^ 9 tanh 



(9) 



- sin^ 6'(a;sin0 - y cos6' + yo cos0)j - iAf sin exp(iMa;). (10) 

This formula describes the oblique solitons generated by the flow of a superfluid past an impenetrable obstacle. It is 
clear from this formula that such solitons can only be generated inside the Mach cone, 



arcsin(l/M) < 9 < arcsin(l/M). 



(11) 



3. If there are several obstacles in the flow of a superfluid, then several dispersive shocks are generated which decay 
far enough from the obstacles into oblique dark solitons. When such space solitons overlap, they interact with each 
other and their behavior in the overlap region is of considerable interest. We have simulated the interaction of oblique 
solitons numerically and the results are shown in Fig. 1. As we see, two pairs of dark solitons are generated, two of 



FIG. 2: Shifts of the obhque sohtons positions as functions of their slope angle. The second soliton has the slope angles: 02—0 
(dashed line); 62 = 0.1 (solid line). The Mach number is equal to M = 5. 



these solitons interact with each other in the region far enough from the obstacles and the end points of the sohtons 
decay into vortices. It is remarkable that the interaction is practically elastic — no new solitons or radiation are visible. 
The only visible result of the interaction is a shift of the solitons positions after their interaction. This behavior is 
typical for the systems described by so-called completely integrable evolution equations (see, e.g. Q). Although there 
is nothing known about complete integrability of the system ([5]), it has well-known limiting cases when it reduces to 
completely integrable equations (see, e.g., 0, 113, first, the limit of shallow solitons, when the system reduces 

to the Korteweg-de Vries (KdV) equation, and, second, the hypersonic limit M ^ 1 when it reduces to ID NLS 
equation. This indicates that the system ([5]) is in a sense "close" to the completely integrable equations and therefore 
it demonstrates similar behavior. In this work, we concentrate only on the study of deep solitons so we shall consider 
the hypersonic limit and derive formulae for the corresponding shifts of the solitons position. 



4. Let us consider the hypersonic limit M ^ 1. As was shown in [T^ [Tl|, in the leading order of the expansion 
with respect to the small parameter 1 /M <C 1 the system ([5]) reduces to 

TiT + {nv)Y — 0, 

^ ^ ^finrr nry\ (12) 

and 

Ui,T + VUi^Y = 0, (13) 

where we have introduced the notation 

M + ui+0{l/M), T=^, Y = y. (14) 

The system (|T2|) is nothing but the hydrodynamic form of the ID NLS equation 

i*T + 5*yy - l^'P* = (15) 

for the variable 



*(r,r) = y/n{Y,T)exp (^i J w(r',T)dy'j 



(16) 



As it is well known, the NLS equation (jlSp is completely integrable, it has exact multi-soliton solutions and interaction 
of two solitons was already studied in the classical paper [1J|. The single soliton solution of the equation is 
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FIG. 3: Left: Cross sections of the density distribution shown in Fig. [T]for different values of x obtained from the numerical 
solution of the 2D NLS equation ([1]). The value of x is indicated on each curve. The collision of two oblique dark solitons 
occurs at a; « 101 and is practically elastic. Right: Cross sections of the correspondent phase distribution. One can see a 
phase jump after the collision of the two solitons . 



parameterized conveniently by the value A of the associated Zakharov-Shabat spectral problem and after returning 
to the a;, y coordinates it takes the form 

1 - 

n{x,y) — l ^ — , (17) 

so that 

X = 9M for|0|<l. (18) 



It is natural that Eq. ([7]) reduces to pT)) in the limit ([ll 

If there are two oblique solitons in the superfluid, then they are characterized by two parameters A1.2 corresponding 
to different angles 6*1^2 — A1.2/M and by two different "initial" coordinates yio, t/20- We suppose that Ai > A2 and 
Vio < 2/20- Then the shifts Ayi of the asymptotic "positions" of the oblique solitons are described by the formulae 

. 1 , (Ai - A2)^ + (^1 + 1^2? . 1 , (Ai - A2)^ + ji^i + 

2vi (Ai - + {vi - 2V2 (Ai - \2Y + [vi - V2y 

where 



X^^Me^, v,^^\-\\, i = l,2. (20) 

These formulae describe the shifts for the case M ^ 1. 

The dependence of the shifts Ayi on 9 [see Eqs. (fT9)) . ((20|) ] for some values of the slope angle of the second sohton 
and i\/ = 5 is shown in Fig. 2. 



6. Now, we compare our analytical predictions with numerical simulations. For large obstacles, many pairs of 
solitons can be generated at different angles past each obstacle [4j . For sake of simplicity, we consider here only small 
obstacles with radius r '-^ 1, thus each obstacle only generates one pair of oblique dark solitons with angles 6 and 
—9. We present in Fig. [3] (left) cross sections of the density distribution shown in Fig. [l]and the correpondent cross 
sections of the phase (right). The collision occurs at x « 101 and is practically elastic, i.e., we do not see any radiation 
loss during and after the solitons interaction. We also observe a phase jump after the collision. 

In order to measure the shifts of the solitons positions after their interaction, we have performed two series of 
numerical simulations. Firstly, we have simulated the 2D flow past one impenetrable obstacle of unitary radius placed 
at (a;io,yio) = (0,-15) and for different values of the Mach number M, namely M = 5, 6, 7, 8 and 10. Then, we 
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have measured the coordinates (a;iA,2/iA) of the minimum of the soHton far enough from the region of interaction, 
where subscript "A" denotes simulation with one single obstacle. 

Secondly, we added another obstacle at the position {x2o, 2/20) = (0, 15), repeated the simulation of the 2D superfluid 
flow and measured the coordinates (a:iB,2/iB) of the minimum of the dark soliton, where "B" denotes simulation with 
two obstacles. Thus, for xia = a^iB, the shift is given by Ayi = ?/ib — J/ia- In all simulations, we have applied a 2D 
finite difference method (Crank-Nicolson method) combined with a split-step method and used the spatial grid sizes 
6x — 5y = 0.2 and the time step 5t = 0.01. 

In Table 1 we compare the numerical results with the analytical predictions for the shifts using different values of 
M. The agreement with the analytical results is satisfactory considering the perturbation of the linear waves [l3| on 
the dark solitons and the computational limitation for the numerical simulations. 



TABLE I: Supersonic flow past two obstacles. 



M 


e 


Aj/i 


Ayi Eq.lISl) 


5 


0.1 


0.5 


0.8 


6 


0.08 


0.6 


0.8 


7 


0.07 


0.6 


0.8 


8 


0.06 


0.6 


0.8 


10 


0.05 


0.5 


0.8 



We further simulate the two-dimensional superfluid flow past a single small obstacle with different sizes and notice 
that the amplitudes and the slopes of these solitons depend on the obstacle's sizes. As it can be seen in Table |lTl 
increasing the size of the obstacle, the amplitude of the soliton increases and its slope decreases. Consequently, one 
can investigate the interaction of oblique dark solitons with different amplitudes and slopes considering two obstacles 
with different sizes. In Fig. |4]we show the superfluid flow past two obstacles, one with radius r = 0.6 and the other 
one with radius r = 1. We see that the collision of two different oblique dark solitons is still practically elastic and 
also the phase jump after the collision. 



TABLE II: Amplitude [A) and slope {6) of dark solitons for obstacles with different radius (r) and M = 5. 



r 


0.6 


0.8 


1.0 


1.2 


1.4 


A 


0.57 


0.69 


0.78 


0.81 


0.83 


e 


0.13 


0.11 


0.1 


0.09 


0.08 
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FIG. 4: Left: Cross sections of the density distribution for different values of x obtained from the numerical solution of the 
2D NLS equation ([l}. The value of x is indicated on each curve. We consider the flow past two different obstacles, one with 
radius r = 1.0 at (xiq, yio) = (0, —12) and the other one with radius r = 0.6 located at (2:20, 3/2o) = (0, 12). The collision of two 
different oblique dark solitons occurs at a; « 103 and is still practically elastic. The solitons angles are 9i =0.1 and 62 = 0.13 
and the shifts of the solitons positions are Ayi « Aj/2 ~ 0.5. Right: Density plot of the corresponding phase. We can see a 
phase jump after the collision of the two solitons. . 



7. Conclusions: We analyzed the collision of two oblique dark solitons numerically and by analytical approximations. 
The observed shifts are consistent in magnitude order with the analytical predictions, considering the perturbation 
of the linear waves and the computational limitation for the numerical simulations. During and after the collision we 
have not observed any radiation loss and phase jumps are analogous to those observed in the ID NLS. We conjecture 
that collisions of oblique solitons in 2D NLS may be a completely integrable process in the asymptotic limit. This 
soliton collision might be experimentally observed in different nonlinear media such as an atomic BEC, photorefractive 
crystals and exciton-polariton condensates. 

We thank A.M Kamchatnov for useful discussions. We also thank funding agencies CAPES, CNPq and FAFESF 
(Brazil). 
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